Investigation of the electromagnetically induced transparency 
in era of cosmological hydrogen recombination 



"'^D. Solovyev, '^•'^V. K. Dubrovich and "^G. Plunien 

" V. A. Fock Institute of Physics, St. Petersburg State University, Petrodvorets, 
Oulianovskaya 1, 198504, St. Petersburg, Russia 
^ Max-Planck-Institut fiir Physik Komplexer Systeme, Nothnitzer Str. 38, 

D-01187, Germany 

- St. Petersburg Branch of Special Astrophysical Observatory, Russian Academy 
of Sciences, 196140, St. Petersburg, Russia 
Nizhny Novgorod State Technical University n. a. R. E Alekseev, LCN, 
GSP-41, N. Novgorod, Minin str., 24, 603950 
Institut fiir Theoretische Physik, Technische Universitat Dresden, 
Mommsenstrasse 13, D-01062 Dresden, Germany 

E-mail: solovyev.d@gmail.com 

Abstract 

Investigation of the cosmic microwave background formation processes 
is one of the most actual problem at present time. In this paper we analyze 
the response of the hydrogen atom to the external photon fields. Field char- 
acteristics are defined via conditions corresponding to the recombination era 
of universe. Approximation of three-level atom is used to describe the "atom 
- fields" interaction. It is found that the phenomena of the electromagneti- 
cally induced transparancy (FIT) takes place in this case. Consideration of 
FIT phenomena makes it necessary to update astrophysical description of 
the processes of the cosmic microwave background formation and, in par- 
ticular, Sobolev escape probability. Additional terms to the optical depth 
entering in the Sobolev escape probability are found to contribute on the 
level about 1%. 



1 Introduction 



In view of recent success in the theoretical description and experimental observa- 
tions of the cosmic microwave background (CMB) a detailed analysis of all the 
processes occurring in cosmological recombination era is required. For a precise 
description of CMB in astrophysical investigations the different photon-emission 
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and absorption, photon-electron scattering and etc. processes should be included. 
The spectral characteristics of the emitted radiation due to different phenomena, 
the polarization state and the cross-section are particularly relevant from the astro- 
physical point of view. In standard calculations for such kind of tasks the quantum 
mechanical approach for the isolated atom is applied. 

However, the consideration of the "atom-fields" interaction becomes impor- 
tant in case of astrophysical experiments with accuracy of about 1% [ID, [|2l and 
with the expectation of even increasing the accuracy up to the level of ~ 0.1%. 
Atom's interaction with an external fields can lead to effects like population in- 
version, steady-state solution, line strengths, cross sections, susceptibility and po- 
larization and etc [|3l-[l8l. Apart from other powerful methods (Green function 
approach, for example) the statistical operator theory can be applied for the study 
of the "atom-field" systems. Application of the density matrix formalism in three- 
level approximation seems to be more simple and appropriate in this case. The 
clear description of the density matrix theory and its applications, such as sponta- 
neous emission, line broadening (power broadening and saturation, collision line 
broadening, Doppler broadening and Voigt profile), can be found, for example, in 

m. 

The radiation transfer theory that is usually applied in the research of cosmic 
microwave background was suggested in [fTO^, [fTTI . In particular, it was estab- 
lished that the 2s ^ Is transition is able to substantially control the dynamics of 
cosmological hydrogen recombination. Moreover, distortions of the order 10~^ 
were predicted [10|. Recently, the radiation transfer theory for the recombination 
era of early universe was intensively re-examined due to precise observations of 
the CMB lfT2ll . lfT3l . In ||T4|| . [15] corrections to the ionization history were found 
that exceed the percent-level. In order to achieve this accuracy the multi-photon 
decays in atoms involving set of states should be included in evaluation of CMB 
formation [[T6l - [|24il . As a rule the two-photon decays of the excited states are con- 
sidered for the evaluation of CMB. Two-photon emission processes were evalu- 
ated accurately in [[25l-[l27|. In our recent works [28J, [29] three- and four-photon 
transitions with separation out of the two-photon links were considered from as- 
trophysical point of view. In [|30l , [1311 we considered one- and two-photon tran- 
sitions in external electric fields. Such modifications should have a strong impact 
to the determination of the key cosmological parameters [1321 . 

Generally, the absorption coefficient calculated per atom is used for the study 
of the radiative transfer in spectral lines. The influence of powerful high-frequency 
electromagnetic radiation on the absorption coefficient in the low-frequency line 
in a three-level A atom was considered in [i33i[ . In this paper we consider the other 
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kind of multi-photon process, namely, the electromagnetic ally induced transparency 
(EIT) phenomena. The nature of EIT phenomena can be examined by evaluation 
of the response of the multilevel system to the presence of an extemal radiation 
field. Electromagnetically induced transparency leads to a significant modifica- 
tion of the absorption profile of the system. Description of EIT phenomena for a 
three-level ladder system interacting with two near-resonant monochromatic fields 
can be found, for example, in [|34l - [|36l . A complete and detailed treatment of 
a three-level ladder atom, absorption and emission spectra, as the transient and 
steady-state response of the S-atom are given in [34J. We study the response of 
the three-level ladder system on the extemal fields originating from photons emit- 
ted during the recombination and evaluate the absorption coefficient, which we 
apply to radiation transfer theory. The physics of an atom interacting with photon 
fields can be understood on the basis of the "interfering-pathways" description, 
which corresponds to multi-photon process defined in term of a power series ex- 
pansion over fields amplitudes (see e.g. Refs. [|34l. [|35l and [f36l ). 

Our paper is organized as follows: In section 2 we shall briefly review essen- 
tials of the density-matrix approach when deriving the density matrices element 
P21 and presenting its series expansion. In section 3 we shall define the absorption 
coefficient and illustrate its application to astrophysics on the basis of the [13J. In 
the subsequent sections we shall provide numerical calculations together with a 
discussion of the results. We conclude the paper with a summary. 

2 Three-level E atom and density matrix 

In this section we employ the three-level ladder scheme for the description of 
the hydrogen atom. We assume that the hydrogen atoms formed during the re- 
combination epoch in earlier universe reach their ground states via emission of 
photons of all the spectral lines corresponding to atomic continuum-bound and 
bound-bound transitions. All the emitted photons generate the (coherent) exter- 
nal field environment, which feats back onto the hydrogen atom. Investigation of 
this self-consistent scenario under the conditions of cosmic expansion is the aim 
of the present paper. We confine our consideration to the spontaneous emission 
rates. CoUisional excitation and ionization can be omitted because at the relevant 
temperatures and densities they are negligible for a three- level hydrogen atom 

m. 

To identify the possible effect we discuss atoms subjected to an extemal field 
which mainly consists of two adjacent spectral lines, namely the Lya and Ha 
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lines, with the initial condition corresponding to full population of the atomic 
ground state. Then we use the standard density matrix formalism which can be 
found, for example, in Ref . (36] . 

The generic three-level ladder system is depicted in Fig. 1. Solving the 
density-matrix equations employing the steady state and the rotating-wave ap- 
proximations yields the following set of equations for the relevant density-matrix 
elements Il36l : 

(P22 - Pll) - ^^P3l] 

721 - iAa 

(P33 - P22) + KP31] 
732 - iAfs 

73i-^(A, + A^) ' ^ ^ 

1 

7rf^(fi^P32-fi/3P23)- 

3 

The levels of system are specified with the following hydrogenic states |1) = |ls), 
1 2) = \2p) and |3) = |3s). The frequencies W21 and ^32 correspond to the transi- 
tions |1) — i- |2) and |2) — |3), respectively. The system is driven by a "probe" 
field with amplitude Ea at frequency Ua and by the "control" field with amplitude 
Ep and frequency up. Aq, = — ijJ2i, A^ = — (X'32 define the corresponding 
de-tunings, Vtp = 2dj,2Ep and Vt^ = 2d2iEa are the Rabi frequencies, which de- 
pend on the atomic dipole-matrix element dij. All these expressions are given in 
atomic units. Neglecting coUisional de-phasing effects the decay rate is given by 
7ij = (Fj + Fj)/2, where Fj is the natural decay rate of the level \i >. 

In the limit of a weak probe field and pu ^ 1, P22 ~ P33 ~ ( full population 
of the ground state of the atom) the solution of the equations ([T]) for the p2i to first 
order in the probe field and to all orders in the control field was found in [|35l . 



z/2[n^ 
P21 = 

P32 — 

P31 = 
P22 = ■ 
P33 = 



4 



Il36l . But the total solution of Eq. ([B for p2i is 



^ ?(S|-^)-(-S| + 7.-.A0(f + ' 

^ = 5(731 -«(Ai + A2)) , 

5 = + 732 - , 

4i 3 

which can be easily transformed in the approximation of weak fields to the ex- 
pression (2) in [|36l . 

The underlying physics of the systems response on external fields can be un- 
derstood considering the power series expansion of the solution Q for the matrix 
element p2i with respect to the variables fi^ and Vtp. Authors [36| derived a se- 
ries that contains terms up to third order in Vt^ and fi^ at zero detuning. For our 
purposes it is important to keep non-zero detuning and, as before, in the approx- 
imation of weak field expression (|2]) will be expanded into a power series with 
respect Vt^ and Vtp. 

The resulting expression for the p2i simplifies significantly in case of exact 
two-photon resonance, i.e. when the frequencies of two external fields coincide 
exactly with the transition frequency uj^i = E3 — Ei. In this case the equality 
Aq, + = holds and the series expansion looks like 



P21 - 



721 - ^Ac 



+ Tr^^^-:r^+ (3) 



731(721 - ^Aa) r2(721 - iAa] 

^»/4 • (r2 + 73i)(72i - - 2731(732 - ^A^) ^ 

(721 - zAc,)2(732 - iAfs) r27|i 

where the dots imply the higher-order terms in fi^, fi^ and product ^q, ■ fi^. The 
series expansion is done under the conditions Vtp/'yij <^ 1 and ^la/lij ^ 1- 

As it was established in [|36l the common pre-factor in (|3]) corresponds to the 
one-photon absorption processes, while the squared terms are associated with the 
two-photon absorption and subsequent emission processes. The products • n| 
represent the "interfering-pathways" terms, see Fig. 2. Thus the matrix element 
P21 describes the multi-photon processes of the coherent "atom-field" interaction. 
For a detailed analysis of Eq. (|3]) we refer to [,36 J . 



5 



3 Absorption coefficient and Sobolev escape proba- 
bifity. 

The theory of radiation transfer for mukilevel atoms utilizing the concept of the 
Sobolev escape probability has been described in [13|. With the method of escape 
probability a simple solution to the radiative transfer problem for all bound-bound 
transitions can be found. The Sobolev escape probability pij (j refers to the upper 
level and i to the lower level of a multilevel atom) is the probability that photons 
associated with this transition will "escape" without being further scattered or 
absorbed. If Pij = 1, the photons produced in the line transition can escape to 
infinity - they do not give rise to distortions of the radiation field. If pij = 0, 
the photons can not escape to infinity; all of them get re-absorbed, and the line is 
optically thick. In general, pij <^ 1 for the Lyman lines and Pij = 1 for all other 
line transitions. The Sobolev escape probability is included in direct astrophysical 
equations (see, for example, Eq. (25) [13]) of radiation transfer. 

Following section 2.3.3 of Ref. [il3il the Sobolev escape probability can be 
presented in the form: 

l-eM-rs) ^ (4) 

where ts is the Sobolev optical depth. The optical depth is a measure of the 
extinction coefficient or absorptivity up to a specific 'depth'. In other words the 
optical depth expresses the quantity of light removed from a beam by scattering 
or absorption during its path through a medium. The Sobolev optical depth can be 
defined as 

rs = ^, (5) 

\v'\ 

where k is the integrated line absorption coefficient and Ajj is the central line 
wavelength. The monochromatic absorption coefficient or opacity is A; = k(j){uij) 
(uij is the frequency for a given line transition and (pii'ij) is the normalized line 
profile), v' is the velocity gradient which is given by the Hubble expansion rate 
Hiz). 

The absorption coefficient depends strongly on the external conditions and re- 
quires the particular consideration for the each case. In presence of an external 
field the opacity can be related to the imaginary part of the density-matrix ele- 
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ments pij as follows: 



(6) 



where is the permittivity of the vacuum and N is the number of atoms. 

Using expression (|3]) for the definition of the imaginary part of p2i, we obtain 



Im{p2i} 



721 

+ 7^1 



[1 + fi^. 



(7) 



together with the dimensionless function 



^^^^^ 



A? 



7ii 



A^ + 7li 



02 

"/3 



O? 



4721731 412721 
(7li732 - 732A^, - 27I17I2 + 2721A2 (37i2 + A/3(3A^ - A,))) 



+ 



16r272i73i(A2+7|,)2(A2 4 
, (7^1 732 - 732A2 



7I2) 

- 272iA,A^) 



(8) 



1672i7|i(A2+72j(A2 



7I2) 



+ 



To get informations about the physical meaning of the function /, i.e. about 
the absorption and subsequent emission processes, one can proceed as earlier for 
the density-matrix element p2i- For the definition of the integrated line absorption 
coefficient from Eqs. (6)-(8) the line profile is separated out. It is assumed that 
the line profile appearing in Eq. (7) corresponds to the monochromatic absorption 
coefficient, see Eq. (31) in [|T3l . Here the function / depends on the fixed param- 
eters Aa and A^, although the common factor represents the Lorentz line profile, 
where Aq, = Ua — ^21- 

Thus the integrated line absorption coefficient can be presented in the form: 



k 



21 



[1 + /(02,02,A„,A^ 



(9) 



with the line profile 0(1^21) = 721/ [{ooa — ^21^ + lii]- accordance with the 
theory described in [[T3l the line profile should be normalized within the interval 
[0, 00] and the coefficient vr arises in Eq. (|9]l. Thus the Sobolev escape probability 
extends to the expression 



P12 



l-exp(-rs[l + /(02,02,A, 



nA,)]) 



/(02,02,A,,A^)] 



(10) 
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where ts can be taken in the form (|5]), Eqs. (39), (40) in [fT3l . In principle, ex- 
pression (flOl ) should be employed in further astrophysical evaluations. We reserve 
this for forthcoming research and restrict ourselves to the consideration of the 
function f{0,'^, i^^, A^, A/j). The function / depends strongly on the parameters 
fi^, fi^, Aa, A/3. The applicability of the power series expansion Eq. (|3]) with 
respect to Qq,, is limited by the field amplitudes and can be found in [35] but 
beyond in ||36l . Estimates for the field amplitudes can be deduced from the CMB 
distribution corresponding to the hydrogen recombination era in earlier universe. 



4 Numerical results and discussion 

While the universe further expanded and cooled down the electrons and protons 
tended towards formation of hydrogen atoms. The temperature at this era is very 
well known from the laboratory physics, T ^ 4500 — 3000 K. After "recom- 
bination", the photons released were able to travel through the universe rela- 
tively undisturbed, and formed the primordial background radiation. However, 
such a photon environment (background) should have influence on the hydrogen 
atom. The field amplitudes for a circular polarized wave can be obtained from the 
(thermal-averaged) spectral energy density 



(11) 



QkgTe — 1 



where c is the speed of light, ks is the Boltzmann constant, h is Planck's constant 
and in further calculations we use Tg = 3000 K. The right-hand side of the equa- 
tion above corresponds to the black-body distribution of the CMB, while left-hand 
side defines the (electrical) energy density. 

In order to avoid de-phasing problem we should choose Az/y ~ F^. Hence for 
the spectral lines 1^21 = i^a (Lya line) and 2/32 = 1^(3 (Ha line) we obtain 

Ea ^ 0.000068802 V/m = 1.33799 • lO^^^a.u., 

^ 52.8636 V/m = 1.02803 • 10^^° a.u. (12) 



The magnitudes (1121) of the field are small; we should compare the Rabi fre- 
quencies with the corresponding level widths. The one-photon transition rates, 
which yield the major contributions to level widths, can be easily evaluated and 
are well known. For the hydrogen atom the dominant one is Ly-a transition rate, 
F2p ~ 10~^ in atomic units and, therefore, the power series (|3} is valid. Moreover, 
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the estimates given in Eq. (fT2l) reveal that we can neglect all terms of the order 
and higher. 

Further we evaluate the function /(f^Q, f^^, Aq, A^) numerically. In case if 
detuning are defined by A^ = and A/3 = + T^g and using Eq. ^ we can 
present the function / in the form: 

f{nl, n% Aa, Ap) ^ -1.30494 ■ lo^^nl - 2.08127 ■ 

+8.67218 • lO^^fi^ + 5.15573 ■ lOV^ + 2.42665 ■ 10^^(]^ (13) 
+3.22916 ■ 10^9^2 _ 1.08979 ■ lO^^fi^^J - 6.18661 ■ lO^^l^^f^J 
-1.63588 ■ lO'^^fi^fiJ + 2.92629 ■ 10^^^]^^]J - 4.88545 ■ 
Taking into account the estimates (fT2l) we receive 

fi„ ^ 1.99343- lO^^^a.M., 
^ 1.92942 ■ 10-^^ a.u., 

f ^ 0.0000322844. (14) 

In case of exact resonances from Eq. ^ and (fT2] ) we get 

A„ = A/3 = , 

/ ^ -0.015814. (15) 

Finally, in case of an exact two-photon resonance, where A^ + A^ = and 
together with Aq, = T2p, the function / takes the value: 

Aa = Tap = -A/3 ~ lO^^a.M. 

/ ^ 0.00952743. (16) 

Thus the magnitude of the function /(l^^, fl'p, A^, A/3) is of about 1% (10"^); 
in case of an exact one-photon resonances of about 1.5%, respectively about 
0.95% if the detuning are no-zero but of opposite sign. In quantum optics this 
effect is well known when the electromagnetically induced transparency is inves- 
tigated for different kind of systems (two-, three- or four level systems with A-, 
V- or H-scheme of levels). The results of calculations for different values for de- 
tuning are compiled in Table 1 and shown graphically in the figures 3 and 4. In 
astrophysical investigations of the CMB formation processes the Lorentz line pro- 
file is used for defining of the absorption coefficient. This line profile represents 
the dominant term in /m|p2i} and is separated out in Eq. (|7]). But in addition 
the function /(fi^, fi^, A^, A^) in Eq. (|7]) has to be included in astrophysical 
evaluation of CMB. 
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5 Conclusions 



The aim of our paper is to study the EIT phenomenon and its influence on the 
CMB formation. The phenomenon of EIT consists in the investigation of the 
systems response on the external field. During the recombination era the hydrogen 
atoms reach their ground state via emission of photons from all the spectral lines. 
After recombination photons were able to travel through universe and formed the 
CMB. Generated by this way the field should affect the primordial atoms. Thus, 
we are confronted with the consideration of the phenomena of electromagnetically 
induced transparancy . 

We employed the quantum optical evaluation of the integrated line absorp- 
tion coefficient that assumes application of photon beams (laser). In astrophysics 
the photon beam diffusion is accounted for via the Sobolev optical depth. Usu- 
ally, the standard line profile is employed for the solution of the recombination 
problem and description of the CMB formation processes and the integrated line 
absorption coefficient is defined via Einstein coefficients in frames of the Sobolev 
approximation. But the determination of the integrated line absorption coefficient 
based on the quantum optical techniques allows one to take into account the in- 
fluence of the external fields on the atom. Accordingly, we derive the absorption 
coefficient from the imaginary part of the density matrix element and we deduce 
the magnitudes of fields from the CMB distribution. 

We are let to the additional function / which depends on extemal conditions. 
The values of function / are listed in Table 1 for different values of detuning. 
The dependence on detuning is depicted in figures 3 and 4. The maximal value 
of the / was found and amounts to about 1.5% in case of exact resonances (when 
both de-tunings are equal to zero). In case of exact two-photon resonance, when 
frequencies of fields are close but differ slightly to the corresponding resonances 
and the total detuning is equal to zero, we obtained / ~ 0.95%. Figures 3 and 
4 also reveal that additional transparency of the medium yields contribution at 
the level about 1%. We expect that the modifications of this magnitude should 
definitely be relevant in determinations of the key cosmological parameters. 

The problem of dephasing appears when defining the magnitude of the field. 
The effect of dephasing leads to the spectral line broadening. To prevent dephas- 
ing phenomena we confined the definition of the field by a narrow strip. In this 
case the width of the corresponding line appears as a natural parameter. In our cal- 
culations we used the following relation Ai/j^ = Ao,(A^). In addition we should 
note the exponential behavior over temperature for the field amplitudes. With in- 
creasing of the temperature Te larger values for the amplitudes could be obtained. 
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see Eq. (fTTI) . Accordingly, the contribution of the EIT effect will become more 
significant for higher temperatures. 
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Figure 1: Scheme of the three-level ladder system under consideration. The three 
levels correspond to the hydrogenic states: |1) — |ls), |2) — )• \2p) and |3) — >■ 
1 3s), respectively. The frequencies cua, coi3 are the frequencies of external fields 
(which correspond to probe and controlled laser fields in [36] ). The fields Ea, 
stimulate transitions Is — 2p and 2p — 3s (Lyman-a and Balmer-a lines). Possible 
detunings and for the amplitudes of the fields are also indicated. 
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Table 1. The numerical results of the function f{Vla, ^p, Aq, for the different 
magnitudes of detunings are presented. In the first column different values for 
/(fifj, Aq, A^) are listed, in the second and third column the detunings and 
field amlitudes are compiled. 



/(fi,,fi^,A„,A;3) 


A„ s-^ 
A^ 

Aui, = |A„|(|A^|)inEq. ^ 


\Ea\ V/m 
\E^\ V/m 


6.48635 ■ 10-^ 


W21 = 6.26826 ■ 10^ 
W32 = 6.31696 ■ 10*^ 


6.8802 ■ 10-5 
5.30686 


0.0000326097 


1^21 = 6.26826 ■ 10^ 
1^21 + 1^32 = 6.33143 ■ 10^ 


6.8802 ■ 10-5 
53.1293 


0.00952743 


1^21 = 6.26826 ■ 10^ 
-W21 = -6.26826 ■ 10^ 


6.8802 ■ 10-5 
52.8636 


0.00952743 


-W21 = -6.26826 ■ 10^ 
1^21 = 6.26826 ■ 10* 


6.8802 ■ 10-5 
52.8636 


0.0000324474 


-W21 = -6.26826 • 10* 
-W21 = -6.26826 • 10* 


6.8802 ■ 10-5 
52.8636 


0.0000326097 


-W21 = -6.26826 • 10* 
-1^21 - W32 = 6.33143 • 10* 


6.8802 ■ 10-5 
53.1293 


4.06852 ■ 10-6 


101^21 = -6.26826 • 10^ 
10(1^21 + W32) = 6.33143 ■ 10^ 


2.1757- 10-^ 
168.01 
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Figure 2: The transition processes occurring in the three-level ladder system cor- 
responding to the different term of Eq. <^ are depicted: part a) of the figure 
corresponds to the one-photon absorption processes (first term equal to the com- 
mon factor); part b) and c) represent the second and third (product) terms and part 
d) describes the "interfering-pathways" (fourth term of |3)). 




Figure 3: The dependence of the function fi^, Aq,, A^) on the detunings 

Aq and A^ is depicted for fixed values of the external fields (fT2] ). i.e. i^p/'jij <^ 1 
and ^la/lij ^ 1- The detunings Aq, and A^ vary within the range [—T2p, ^2p]- 
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Figure 4: The section through the 2-dimensional plot of function 
/(fi^, Aa, A^) (see [3]) at A^ = is plotted. The remaining detuning 
parameter A^ varies within the range [— 3r2p, 3r2p]. 
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